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ABSTRACT 



Asteroseismology with the space-borne missions CoRoT and Kepler provides a powerful mean of testing the modeling 
of transport processes in stars. Rotational splittings are currently measured for a large number of red giant stars and 
can provide stringent constraints on the rotation profiles. 

The aim of this paper is to obtain a theoretical framework for understanding the properties of the observed rotational 
splittings of red giant stars with slowly rotating cores. This allows us to establish appropriate seismic diagnostics for 
rotation of these evolved stars. 

Rotational splittings for stochastically excited dipolar modes are computed adopting a first-order perturbative approach 
for two 1.3A/o benchmark models assuming slowly rotating cores. 

For red giant stars with slowly rotating cores, we show that the variation of the rotational splittings of ^ = 1 modes 
with frequency depends only on the large frequency separation, the g-mode period spacing, and the ratio of the average 
envelope to core rotation rates {TZ). This leds us to propose a way to infer directly TZ from the observations. This 
method is validated using the Kepler red giant star KIC 5356201. Finally, we provide a theoretical support for the use 
of a Lorentzian profile to measure the observed splittings for red giant stars. 

Key words. star:evolution- stars: interiors- stars: rotation- stars; oscillations 
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1. Introduction 

Stellar rotation plays an important role on the structure 
and evolution of stars. The problem of transport of angular 
momentum inside stars is not yet fully understood, how- 
ever. Several mechanisms seem to be active although they 
are either only approximately described or not modeled at 
all. In order to study the internal transport and evolution 
of angular momentum with time, one needs observational 
constraints on physical quantities affected by such trans- 
port processes. In particular, the knowledge of the internal 
rotation profiles and their evolution with time is crucial. An 
efficient way is to obtain seismic information on the internal 
rotation profile of stars. 

The ultra-high precision photometr y (UHP) asteroseis - 
mic space missions such as CoR oT (jBaglin et all I2006D 
and Kepler (jBorucki et al.l |2010() offer such an opportu- 
nity. Because cool low-mass stars have a convective en- 
velope, oscillations can be stochastically excited as in the 
solar case. Red giant stars are of particular interest here. 
Stochast ically excited non-radi al modes were detected with 
CoRoT (|De Ridder et al.ll2009D . Analyses of these data re- 
vealed the oscillation propertie s of a large number of these 
stars (e.g., 



isciiiation propertie s or a large numper oi tnese 
Hckkcr e t al. l l2009t iMosser et al.l[20Tnl l2011al: 



iBedding et al . 2011). Their frequency spectra show similar- 
ities as well as differences with the solar case (jMosser et al.l 
l2011bt) . Several works have investigated the properties of 



the non-radial oscillation mode s of red giant stars (e.g. , 
Osakil [T975I: iDzienibowskil [19771 iDziembowski et al] 120011: 



Dupret et al.l l2009t iMontalban et al.l |2010[) . The internal 
properties of a red giant star, a dense core and a diffuse 
envelope, give rise to both gravity-type oscillations in the 
central region and acoustic-type oscillations in the enve- 
lo pe. These o s cillati ons were first called mixed modes by 
Dzi embowskil (|197lD for a Cepheid model and IScuflairel 
(1974y for condensed models representative of a red giant 
structure. This has already been pr oved useful to probe 
the structure of th ese stars (e.g., IBedding et al.l 120111 : 
IMosser et aIll2011aD . 

More recently, measurements of rotational splittings of 
red giant stars ha ve been obtained by us i ng the Kepler ob- 



servations (e . g., Deheuvels et al.l l2012t iBeck et all 120121 : 



IMosser et all 1201 2bD . These splittings provide the first 
direct insight into the rotation profiles of the inner- 
most layers of st ars ([Deheuvels et al.ll2012|) . Observations 
(Mosser et al. 20 12al) reveal that a sub-sample of oscillating 
red giants have complex frequency spectra. Their central 
layers are likely to rotate quite fast. A correct rotation rate 
ought then be deduced from direct calculations of frequen- 
cies in a non-pertubative approac h (see the case for poly- 
tropes and acoustic modes in , Reese et al.l l2006l: Lignieres 



2011 



2010 



lOuazzani et al.|[2012l and for g modes in [Ballot et al 
201 it) . This will be necessary in order to decipher 
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the complex spectra for fast rotation. On the other hand, 
the rest of the sample shows simple power spectra that can 
be easily identified without ambiguity. These ide ntifications 
then lead to quite small rotational splittings (jBeck et alJ 
I2OI2I : iMosser et al.ll2012al ). thus to slowly rotating cores. 

On the theoretical side, the central rotations pre- 
dicted theoretically by standard models currently includ- 
ing rotationally induce d mixing are far too fast (e.g., 
lEggenberger et aLll2012[ ) and cannot account for these ob- 
served slow core rotation in red giants. 

Motivated by these recent results, we have started a 
series of studies dedicated to an understand ing of the rota- 
tional properties of oscillating evolved stars. iMaraues et alJ 
(j2012[) (hereafter Paper I) computed the rotation profiles 
and their evolution with time from the pre main-sequence 
to the red giant branch (RGB) for 1.3 Mq stellar mod- 
els with an evolutionary code where rotationally induced 
mixing had been implemented. The authors then followed 
the evolution of the rotational splittings calculated to first- 
order approximation in the rotation rate using the rotation 
profiles predicted by their evolutionary computations. The 
authors found that for low-mass stars, the first-order ap- 
proximation is valid for most stochastically excited oscil- 
lation modes for all evolutionary stages except the RGB. 
For red giants, the authors showed that there is indeed 
room enough in the uncertainties of the description of trans- 
port of angular momentum to decrease sig nificantly the 
core rotation of red giant models (see also iMevnet et al.l 
l2012f) . However, they still predict a core rotation that is still 
too large compared with the observed rotational splittings, 
meaning that some additional processes are operating in the 
star to slow down its core. To test such possible additional 
mechanisms, investigations of the specific properties of the 
rotational splittings of red giants with slow core rotation 
must be carried out, essentially because of the particular 
dual nature of the excited red giant modes. 

We thus defer the study of red giants with fast rotating 
cores to the third paper of this series. We focus here on 
red giants with slowly rotating cores. We present a study 
of the properties of linear rotational splittings, computed 
for adiabatic oscillation modes, in the frequency range of 
observed modes in red giant stars. We assume arbitrarily 
slow rotation profiles and then study essentially the impact 
of the properties of the modes on the rotational splittings. 
Our goal is to provide a theoretical framework for the inter- 
pretation of the observed rotational splittings of red giants 
in terms of core and envelope rotation. 

The paper is organized as follows; in Sect. [21 we de- 
scribe the properties of our models. In order to interpret 
the rotational splittings, we need first to identify the phys- 
ical nature of the excited modes. This is done in Sect. [3] 
where we distinguish two classes of mixed modes. In Sect, 
m we compute theoretical rotational splittings using adi- 
abatic eigenfrequencies and eigenfunctions obtained for £ 
= 1 modes. The studied frequency range corresponds to 
expected stochastically excited modes. We then show that 
the rotational splittings exhibit repetitive patterns with fre- 
quency that can be folded into a single pattern as they 
carry nearly the same information. We then consider in de- 
tail such pattern for one of our models and investigate the 
information provided by the corresponding rotational split- 
tings. In Sect. [5j we establish an approximate formulation 
for the linear rotational splittings in terms of the core con- 
tribution to mode inertia. The approximate expression fits 
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Fig. 1. Hertzsprung-Riissell (HR) diagram showing a 
I.3M0 evolutionary track. The red crosses indicate the se- 
lected models as described in Sect. [Hand in Table. [1] 



Table 1 

models. 



Fundamental parameters of the 1.3 Mq stellar 



model 


Ml 


M2 


R/Rq 


3.58 


6.91 


logTeff (K) 


3.694 


3.758 


log L/Lq 
Pc (10^ kg/m^) 


0.825 


1.345 


7.06 


14.05 


P./P (10') 


5.12 


75.85 


rcz{mcz) 


0.25 (0.33) 


(0.078) (0.197) 




0.01 (0.13) 


(4.5 10"^) (0.168) 


Al/ (mHz) 


23.12 


8.47 


An (s) 


96.64 


71.14 




342.6 


83.5 


^max 


14 


9 



well the numerical ones and allows us their interpretation. 
We also give a procedure that enables ones to derive the 
mean rotation of the core and of the envelope. Proceeding 
further in the approximations in Sect(5J we find an approxi- 
mate formulation of the core contribution to mode inertia - 
thus to rotational splittings ~ that depends only on observ- 
able quantities aside rotation. We validate our theoretical 
results on the Kepler observations of the star KIC 5356201 
(Beck et al. 2012). Finally, we end with some conclusions 
in SectE 



2. Structures of stellar models and rotation profiles 

We consider 1.3Mq models of red gi ants as this mass i s 
typical for observed red giant stars (jMosser et al.l l20l"ol ). 
These models (structure and rotation profile) have been 
evolved from PMS as described in Paper I. 



2.1. Selected models 

The first model (Ml) lies at the base of the red giant branch 
(Fig.[T]). Its outer convective zone occupies 67% in mass of 
the outer envelope and is located at Am/M = 0.2 above 
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Fig. 2. Logarithm of normalized squared Brunt- Vaisala iV^ 
(solid curve) and {i = 1) Lamb frequencies Sf (dash- 
dotted curve) as a function of the normalized mass m/M 
for model Ml (logj^Q cr^ax = 2.60) (black curves) and M2 
(logio ^max — 2.32) (blue curves). The thin horizontal solid 
lines indicate the frequency range for expected excited os- 
cillation modes: black for model Ml and blue for M2. 

the H-shcU burning. The second model (M2) is further up 
on the ascending branch (FigU]). Its outer convective zone 
occupies 80 % of the outer envelope and is at Am/M — 
0.015 above the H-shell burning. Table [T] lists the stellar 
radius in solar unit, the effective temperature, luminosity 
in solar unit, the central density pc, the ratio pd P with p = 
3M/(47ri?^), the radius at the base of the outer convective 
zone, the radius of maximum nuclear production rate, both 
normalized to the stellar radius, the seismic quantities: the 
large separation Ai/, the g mode period spacing All (in s), 
the frequency at maximum power spectrum intensity Vmax, 

the radial order rimax at t'max- 

The propagation diagrams (|Osakilll975f ) for models Ml 
and M2 are displayed in Fig. [51 Only the central part is 
shown as it plays a major role in determining the properties 
of the oscillation modes of red giants (see Sect. [3]). The 
Brunt- Vaisala {N) and Lamb {Si) frequencies take their 
usual definition, i.e., 

y2 ^ 3 ( ^ ^^'^P dlnp \ 
r yFi dlnr dlnry ' 

Sf^ = , (2) 

where A = t{i -f 1), £ is the angular degree of the mode, 
is the squared sound speed, the gravity g — Gm/r^ , and 
p, /9, Fi have their usual meaning. Unless otherwise stated, 
we consider all squared frequencies iV^ , Sf to be normalized 
to GM/R^ and cr^ = uj'^/{GM/R^) (where uj is the mode 
pulsation in rad/s). 

Model M2 is more evolved and more centrally condensed 
than Ml. Accordingly, the inner maximum of the Brunt- 
Vaisala frequency is hig her than for the younger model Ml 
(jMontalban et al1l2010[ ) and N decreases more sharply. The 



Fig. 3. Rotation profiles as a function of the normalized 
mass. The surface values are f2surf/(27r) = 0.154, 0.042 /xHz 
respectively. Original 17 profiles stemmed from the evolu- 
tion of the stellar models (red solid curve for model Ml 
and blue solid curve for model M2). The dot-dashed lines 
represent the rotation profiles that have been arbitrarily 
decreased by assuming ^suA^/^surif'^ with / = 3.85 for 
model Ml (red curve) and / = 2.8 for model M2 (blue 
curve) . 

range of frequency corresponding to the expected excited 
modes are indicated by horizontal lines. They show that 
the evanescent region between the inner g resonant cavity 
and the upper p resonant cavity {i.e., between the Brunt- 
Vaisala and Lamb frequencies) is quite narrow for these 
modes. 

2.2. Rotation profiles 

The original rotation profiles obtained for models Ml and 
M2 from evolution are displayed in Fig. [31 The central re- 
gions of the more evolved model M2 are rotating much 
faster than those of Ml (amounting to 181.5 and 260 pRz 
respectively). With such high rotation rates, we checked 
first that the centrifugal acceleration is negligible compared 
to local gravity in our models (that is, with a relative mag- 
nitude of 10~^ at most). We can therefore still safely as- 
sume a spherically symmetric equilibrium state. However, 
the linear approximation is no longer valid due to Coriolis 
effects when 2fl/uj is close to unity. This happens in the 
central regions of the models for the modes we consider. As 
we study here slowly rotating red giants, we assumed that 
some slowing down process has and/or is occurring in real 
stars that are not yet included in our models. We there- 
fore artificially decrease the rotation rates stemmed from 
our evolutionary models so that the first-order approxima- 
tion for the rotational splittings is valid. The shape of the 
rotation profile due to slowing down is unknown, we there- 
fore keep it unchanged and investigate whether this can be 
confirmed or not by confrontation with observations. 

We therefore decrease the central rotation rates by 
rescaling the rotation profile stemmed from our evolution- 
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ary models, decreasing the central values but keeping the 
same envelope rotation rate. The decrease is set so that 
the value of the computed rotational splittings agree with 
observed ones. For instance, for model Ml, the initial cen- 
tral rotation, 181.5 /iHz, is decreased down to 0.967 /iHz. 
The rotation profiles we use to compute the linear rota- 
tional splittings in the following sections are then shown in 
Fig. 13 

3. Classification of mixed modes 

The rotational splittings depend not only on the rotation 
profile but also on the eigenfunctions - that is on the phys- 
ical nature - of the excited modes. Stellar models of red 
giant stars exhibit an outer convective region that is able t o 
efficiently drive modes stochastically (jSamadi et al.ll20O) . 
Thus, we consider the frequency ranges for our models that 
span an interval of a few radial orders below and above 
nmax (the radial order corresponding to the frequency at 
maximum power). We estimate rimax as rimax = t'max/Az^, 
where the frequency at maximum power spectrum J^max and 
the mean large separation Av are given by the usual scaling 
relations (e.g., iKieldsen fc Beddi"n3ll995[ ) 



^max,0 



M 



Av 



1/2 



R 
R 

Rq 



off 



cff.O 



-1/2 



-3/2 



(3) 
(4) 



with Toff,© 

UeS t-max,© ■ ^ ^ 

calibrated by'Mosser et al.' ^2013' 



Osakil 



5777 K and the reference asymptotic val- 
3106 /xHz, an d Avp , = 138.8 /xHz defined and 
Indeed, it has been con- 
jecture d by Brown et al. (1991), then shown observation- 
ally bv lBedding fc Kieldseni (j2003[) . that fmax predicts well 
the location of the excited frequency range of stochastically 
excited mod es. A physical justifica tion has been recently 
proposed bv lBelkacem et all (|201l[ ). The adiabatic oscilla- 
tion frequencies and eigenfunction s for i = 1 modes are 
comp uted with the ADIPLS code (jChristensen-Dalsgaardl 
I2008D . 

'e^,^ 



As expected from previous works , 

1 1 - ■ : ■■■ '1977': 'Dzicmbow ski et al]l200lt iDupret et al.l 
2009 : Mon talban et a l. 2010), our red giant models show 
many mixed £ = 1 modes. For model Ml, for instance, 
the frequency range of excited modes is expected to be 
\ogn^{a^) ^ 2 — 2.8. For modes that are expected to be 
excited around v^ax (given in Table the corresponding 
frequencies are much smaller than both the Brunt- Vaisala 
frequency and the Lamb frequency in the central regions of 
the models (see Fig. Therefore, modes with frequencies 
close to Vmax are trapped in the g resonant cavity, located 
in the radiative central region. Moreover, cr^^ax ^ N'^^Sj 
in most of the envelope so that the modes are also trapped 
in the p resonant cavity in the convective envelope (not 
shown). These cavities boundaries are delimited by turn- 
ing points, i.e. the radii satisfying cr^ = N'^{r) (g cavity 
bounded by xi and X2), or = Sf{r) (p cavity bounded 
by X3 and X4). The location of these turning points for a 
list of selected modes are given in Sect. 4-1, Tabled 

In order to interpret correctly the observed rotational 
splittings, we need to define a measure of the g nature 
of the modes with respect to their p part and classify 
them accordingly. Convenient quantities for that purpose 



^ -4 
S. 

^ -5 
£ -6 

O 

" -7 

-8 
0.5 
_ 0.4 
^ 0.3 
B 0.2 

-t-J 

I 0.1 




n ' ' ' r 

model Ml 
0„„ 
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0.88 /j,Hz 
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14 



— 1 1 [— 



r » * ■» 



10 



12 



14 



Fig. 4. Top: Integrated kernel (/? /) with /3 and / respec- 
tively given by Eq. [T2] and Eq. [5] (blue dotted line, dots 
correspond to modes) and mode inertia / (red dotted line, 
dots correspond to modes) in decimal logarithm for model 
Ml. Black crosses represent the inertia of radial modes. 
Bottom: Rotational splittings (EgfTTj). in /iHz, ioi £ — 1 
modes as a function of the normalized frequency v j Av for 
Ml (black dots connected by a dotted line). Dotted blue 
and magenta curves and blue crosses represent various con- 
tributions and approximations to the rotational splittings 
(see text SectO. The value (l/2)(^2)corc/27r is represented 
by the horizontal red dotted line. 



are the kinetic energy ( Unno et al.l[T989[) or the mode iner- 
tia ([Dziembowski et aLl200lD . defined as 



(5) 



with A — £{£ + 1), and x = r/R is the normalized radius. 
The oscillation quantities entering the above definition are 
the fluid vertical and horizontal displacement eigenfunc- 
tions ^r,Ch respectively. For later convenience, we will use 
the following variables 



Zl 



1/2 



^ ] ^3/2 §1 



2:2 = \/a 



3p 

J) ^ R 

2P] 2-3/2 |h 

pj R 



(6) 
(7) 



Using Eqs. ([5]) and ([7]) into (O leads to (dropping the sub- 
scripts n, £) 

\zl + zi)'^, (8) 

where we have defined the dimcnsionless mode inertia / as 
I = I/{MR^). 

The variation oi £ = 1 mode inertia with frequency is 
shown in Figs. 2] and [5] for models Ml and M2. The number 
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27T 
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10 
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1 



Fig. 5. Same as Fig. g] but for model M2. 
of nodes in the g cavity is given by 




(|Unno et al.lll989l and references therein). 

Because A^^ increases with evolution (see Fig. [3] for Ml 
and M2), the number of nodes, hence of modes, in a given 
frequency range increases with the age of the star. This 
can be seen by comparing Fig. 2] for model Ml and Fig. [5] 
for model M2. Mixed modes with a g-dominated character 
have their inertia larger than modes that have their inertia 
shared between the g and p cavities . These results are in 
agreement with previous works fe.g-. lDupret et al.ll2009t ). 

We now define for convenience a measure of the g nature 
of the mode with the ratio of mode inertia in the g cavity 
over the total mode inertia (jPeheuvels et al.ll2012f) . i.e., 

For a mode with a p-dominated nature, the inertia is con- 
centrated in the p-cavity and is small due to the low density 
in the envelope ~ 0) . Its frequency is mainly determined 
by the p-cavity and is therefore half a large separation away 
from the frequencies of consecutive radial modes. 

For a mixed mode where the g character dominates (C ~ 
1 ), the inertia is large due to the high density of the central 
regions where the mode has its maximum amplitude. The 
frequency of this £ = 1 mixed mode differs significantly 
from that of a £ = 1 mode that would be a pure p mode. 
Its frequency is then closer to that of the closest radial 
mode. We re f er later to such a mode as a g-m mode as in 
iMosser et all (|2012a[) . 

When the contributions to mode inertia from the core 
and the envelope are nearly equal, the mode frequencies are 
less affected by the core and therefore remain closer to the 
frequencies of pure i = 1 p modes. We will refer to such 



Fig. 6. Folded rotational splittings according to {i^ — 
va)/I\v for Top: model Ml (blue 9 < v / I\v < 15; ma- 
genta ly/Av > 15; green v/Av < 9) and Bottom: for model 
M2 (blue 8 < i^/Aiy < 14; magenta i^/Aiy > 14; green 
vj Av < 8). The splittings are normalized to their maxi- 
mum values, which are written in the panels. 

mixed modes as p-m modes (|Mosser et al.ll2012ah . These 
modes take intermediate values of C. 

4. Linear rotational splittings 

For slow rotation, a first-order perturbation theory pro- 
vides the following expression for the rotational splittings 
()Ledouxlll951l : [Christensen-Dalsgaard fc Berthomieulll99lD 

^" = 1! (11) 

where fl is the angular rotation velocity and the rotational 
kernel K takes the form 

+ (12) 

normalized to the mode inertia /. For later purpose, we also 
define ^ 

13= I K{x) dx. (13) 

We compute the rotational splittings according to 
Eq. ([TT|) for £ = 1 modes and will refer to them as numer- 
ical rotational splittings later on. The variations of the ro- 
tational splittings with the normalized frequency v j Av are 
shown for model Ml an M2 in Figs. 2] and [S] Also displayed 
are the variations of both the inertia / and the product /3 / 
(with /3 and / respectively given by Eq.fT^and Eq.[5]). They 
show that the oscillatory variations of the splittings with 
frequency are dominated by the nature of the modes as they 
closely follow those of mode inertia and of the integrated 
kernel /3 I (see Sect.Oj). 
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It may be useful to the reader to note that models that 
are located at the same location in the HR diagram (mod- 
els with the same mass but one with an overshoot of 0.1 
pressure scale height and no rotational induced mixing and 
one without overshoot but with rotational induced mixing) 
have similar structures. As a consequence, the rotational 
splitting of any given mode in the range of interest here is 
the same for both models when computed with the eigen- 
functions of the respective models and the same rotation 
profile. 



4.1. Periodic patterns for rotational splittings 

It is clear from Fig. S] and [5] that the variation of the 
mode inertia and rotational splitting with frequency shows 
a quasi-periodic pattern, with a period close to the large 
separation. Each pattern around a radial mode carries 
roughly the same information about the rotation, at least 
in the range o f obser ved modes for red giants as found by 
osser et all (|2012al) . Note, however, that at higher fre- 
quencies, mixed modes become closer to pure p modes and 
probe the rotation of the surface layers, while at lower fre- 
quency, the nature of g-m modes is much closer to that of 
pure g modes and probe better the inner rotation. 

The existence of repeated patterns leads us to gather in- 
formation by folding the power spectrum according to the 
large separation Av. The rotational splittings are plotted, 
in Fig. [HI as a function of {ly — vq)/ Av where z/q is the clos- 
est radial mode to the I — 1 mode of frequency v. The 
folding {v — i'o)/Av corresponds t o the asymptotic fold - 
ing v/Av — {up + e) as suggested bv lMosser et al.l (|2012bD . 
where Up is the radial order of the mode and e r epresent 
some departure from the asymptotic description (jTassoull 
119801: iMosser et aLll2012aD . This folding is shown for our 
two models in Fig. El We see that all patterns do nearly 
superpose for each model. Here e ~ 1.3 — 1.4 for model Ml 
and e ^ for model M2. A clear correlation exists between 
the magnitude of the splitting and the g nature of the mixed 
mode. The rotation splittings are large for g-m modes that 
lie on the left and right sides of a minimum occupied by 
the p-m modes that have the smallest rotation spl ittings. 
This is in agreement with recent observations ( Moss er et al.l 
l2012bl) . In the low- frequency domain, the modes are, for 
most of them, g-m modes. At higher frequency, they can be 
either g- or p-m modes (Fig. [5]) . This is clearly seen when 
one looks at restricted domains of frequency (corresponding 
to different color codes in Fig IB]). 

We also note that the ratio of the minimum splittings 
'^^'min to the maximum ones 5fmax (~ 0.5 for for model Ml 
and M2) is only slightly sensitive to rotation. Actually, it 
depends only on the ratio of the envelope to core rotation 
which is small. This is explained in SectlHl 



4.2. Properties of a periodic pattern 

To interpret the behavior of the splittings, we consider a 
single pattern around a given i — mode with radial-order 
Up = 10 (Table[2]) for model Ml. There are six modes in the 
pattern labelled i^i to i^^ from the smallest to the largest 
frequency. 



Table 2. Model Ml frequencies of the selected pattern: 
the frequency v in /iHz, cr is the normalized frequency 
2'Ki>{GM/R^)~'-^/^ , I is the mode inertia arbitrarily normal- 
ized to mode inertia of the radial mode n = 8; — Icorc/I 



£^0 


rip 


u 


a' 


vjAv 


C 


7(10-") 


10 


261.33 


233.33 


11.31 


6.09 10"" 


0.375 


e = 1 


freq. number 


V 




vjl^v 


c 


7 (10-") 


#1 




212.88 


10.80 


0.422 


0.867 


#2 


252.93 


218.58 


10.94 


0.799 


2.333 


#3 


258.52 


228.34 


11.18 


0.937 


7.512 


#4 


264.70 


239.37 


11.45 


0.928 


5.759 


#5 


270.52 


250.02 


11.70 


0.657 


0.922 


#6 


273.65 


255.84 


11.84 


0.498 


0.578 



Table 3. Normalized radius of the turning points for se- 
lected modes for model Ml frequencies: the g-cavity is 
delimitated by [a;i,X2] and the p cavities by [a;3,a;4]. In 
the other regions the modes are evanescent. We find x\ = 
1.3 10-"* and X4 = 0.9999 for ah listed modes 





X2 


X3 


#1 


0.0732 


0.119 


#2 


0.0725 


0.118 


#3 


0.0715 


0.116 


#4 


0.0705 


0.113 


#5 


0.0693 


0.111 


#6 


0.0689 


0.110 



4.2.1. Mode inertia / and rotational splittings 

The rotational splittings for vi to fg are displayed in Fig|7l 
The relative contributions Ii/I related to zf, and I2/I re- 
lated to z|, defined as 

are displayed in Fig IT] 

The contribution of the Z1Z2 terms to the rotational 
splittings (Eg fTTl) is negligible in front of zf in the envelope 
and in front of z| in the core. Thus, the following conclu- 
sions apply for both the mode inertia and the rotational 
splittings. 

The Ii contribution coincides to a good approximation 
with the envelope contribution Icnv, while the I2 contribu- 
tion coincides to a good approximation with the core con- 
tribution. The core contribution to inertia, /core, and the 
rotational kernel, Kq^cotc = /co^e ^(■^) ^i^) da;, are com- 
puted between the lower and upper turning points of the 
g cavity for each mode. The dense core contributes heavily 
to the total mode inertia / and rotational splittings. The 
values of the ratio C = Imic/I are given in Table [51 

— The modes with frequencies and 1^4 are g-m modes. 
For these modes, the contribution of the horizontal dis- 
placement Z2 (essentially in the core) dominates their 
inertia {I2/I ^ Icoic/I ~ 1) while the contribution due 
to the radial displacement zi (largest in the envelope), 
/i//, is negligible. 

The envelope above the core where the mode is either 
evanescent or of acoustic type contributes for almost 
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Fig. 7. Top: relative inertia contributions Ii/I (green dots 
connected dashed line dots correspond to modes) and I2/I 
(Eq fH|) (black), relative contribution C, = I core /I (magenta 
dashed line, hexagons correspond to modes). Bottom: con- 
tribution to the rotational splittings for model Ml due to 
z1 only (green dashed line dots correspond to modes), con- 
tribution due to z| only (black dashed lines and hexagons) , 
contribution due to Z1Z2 only (red-dashed line, dots corre- 
spond to modes) ; contribution to 5vni due to the core only 
(magenta filled hexagons) 

nothing to the total mode inertia Icm/I = 1 — C ^ C 
for such modes. The contribution of the envelope to Sv 
is negligible because the rotation is quite small there. 

— The other modes are p-m modes. The envelope above 
the core where the mode is either evanescent or of acous- 
tic type contributes for half for these modes (/i ~ /2 ~ 
J/2). They share their inertia almost equally in the core 
and in the envelope with an almost equal contribution 
of the horizontal and radial displacements. 



4.2.2. Integrated kernels 

The modes are trapped in the same g cavity, thus they 
probe the same inner part of the rotational profile of the 
model. This can be seen with the behavior of the integrated 
rotational kernels that are shown in Fig. [5] (left panel). For 
all modes of the pattern, the integrated rotational kernels 
increase rapidly from the center. 

— For the g-m modes, the integrated kernels saturate quite 
rapidly to the maximum value, indicating that the en- 
velope does not significantly contribute. This is the case 
for 1/4 in our studied pattern. Note that the saturated 
value is obtained long before the upper turning point of 
the g cavity is reached. This is due to the fact that the 
amplitudes of the displacement become rapidly small in 
the upper part of the g-cavity. 



Fig. 8. Integrated kernels K{x') dx' (left panel) and 
/p^ f2(x) K{x') dx' (right panel) (with K{x) given by 
Eg. 1121) as a function of x = r/i? for some modes of the se- 
lected pattern for model Ml. The lines correspond to mode 
# 1 (blue); # 2 (black); # 4 (red); # 5 (green) and # 6 
(magenta) of Tab [5] respectively. 



The p-m modes such as vi and z^g have their inte- 
grated kernels that saturate less rapidly and to a much 
smaller value than the g-m modes. The corresponding 
integrated contributions to the splittings are shown in 
Fig. [3 and [8] (right panel). As they include the rotation 
profile, they increase faster than the integrated rota- 
tional kernel. 



5. Seismic constraint on rotation 

To have more physical insight into the information carried 
by rotational splittings it is worthwhile to emphasize their 
dependences. To this end, we decompose the rotational 
splittings into two contributions. Given that the inner (ra- 
diative) and outer (convective) cavities have well distinct 
properties, we write: 



where 

/^corc = 
Penv = 

and 

(^) core 



/ n 



K{x) dx , 



2^ 



K{x) dx = 13 - /3core , 



^ J^^-n{x')Kix') dx' 
K{x') dx' 



(15) 

(16) 
(17) 

(18) 
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Jl^^^njx') K{x') dx' 
K{x') Ax' 



(19) 



where K{x) is defined in Eq. (IT^ and (ri)corc is the angular 
rotation velocity averaged over the central layers enclosed 
within a radius Xcorc = fcmc/R corresponding to the g res- 
onant cavity. The radius Xcoro is taken as the radius of the 
upper turning point of the g cavity X2 (see Table [3] for one 
pattern). (i7)cnv is defined as the angular rotation velocity 
averaged over the layers above Xcoro- 

For the modes we consider in this paper, the evanescent 
regions above the g resonant cavity and above the p cavity 
are narrow. Then, the envelope essentially corresponds to 
the p resonant cavity. The boundary radii are then taken 
as ^3 and X4 (see Table [31). Strictly speaking, the above 
quantities depend on the mode (actually on its frequency) 
because they are defined in the propagative cavities char- 
acteristics of each mode. However, because of the sharp de- 
crease of the Brunt- Vaisala frequency at the edge of the H 
shell burning region, the upper turning point of the g cavity 
is roughly the same for all modes (see Table |3]). The modes 
have a significant amplitude in the g cavity in a region that 
is quite narrower than the extent of the g cavity and that is 
independent of the mode (in the frequency regime we con- 
sider). This can be seen in Fig. lA.ll where the horizontal 
and radial displacement eigenfunctions are shown for two 
typical modes of the studied pattern with frequencies listed 
in Table[21 We can then consider that Xcoi-c and the quantity 
(^^)core are nearly the same for all considered modes. 

The radius of the lower turning-point, X3, for the p- 
cavity of p-m i = 1 modes decreases slightly with increas- 
ing frequency of the mode while the upper turning point 
remains approximately the same (see Table. [3]). Thus, mode 
inertia remains dependent on the mode. On the other hand, 
the mean envelope rotation is nearly that given by the uni- 
form rotation of the convective zones and remains the same 
for all modes. 

(17)e„v « flcZ ■ (20) 

In Fig. m the rotational splittings obtained with two 
approximations for (f2)core are compared with the exact 
calculation Eq. (fTT|) . The core contribution to the rota- 
tional splittings (ri)corc/5corc (Eqs. (dni) and (US])) is dis- 
played with the magenta dots connected with the magenta 
dotted line. The blue crosses dotted line represents the core 
contribution to the rotational splittings (Eq. ([TS])') with 
(r2)core computed with only the horizontal eigenfunction 22- 
The two variations (magenta dots and blue crosses) can- 
not be distinguished. The blue dotted line represents the 
contribution from the envelope to the rotational splittings 
(r2)env (/3 — /3core)- This Contribution is much smaller for 
model M2 than for model Ml. 

For the g-m modes (i.e. for 1^ ~ 1), /3corc ~ 1/2. The 
maximum rotational splitting is then given by 



2^ 



(21) 



The maximum values of the computed rotational splittings, 
^t'max, reach 0.44 and 0.430 //Hz for Ml and M2 models re- 
spectively. These values give rise to a mean core rotation 
(r2/27r)^Qj.^ — 0.88 fiRz and 0.86 fiRz, which are close to but 
smaller than the central rotation for models Ml and M2 re- 
spectively (resp. 0.97 and 0.95 /xHz). This is explained by 
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Fig. 9. The rotational splitting normalized to its maximum 
value as a function of ( (black dots) for model Ml (top) and 
model M2 (bottom) . The red crosses represent the approx- 
imation (5i^/(5i^max = C as a function of C- The blue crosses 
represent the approximation Sv/Svaia,^ — — 2TZ) + 2TZ as 
a function of ^ with R — 0.1785 and using numerical values 
ofC 



the fact that the rotation decreases sharply in the central 
regions of our models and that the major contribution to 
the rotational splitting is slightly shifted off the center as 
are the maximum amplitudes of the horizontal displace- 
ment eigenfunctions (Fig. lA.ll) and the rotational kernels. 

Moreover, for ^ = 1 modes, one obtains /Score — C/2 , 
/3 = 1 — C/2, and /3cnv = l~Ci to a very good approximation 
(Eq. (|A.5p in Appendix IA.2[) . Then the linear rotational 
splitting Eq. (fTS)) is easily rewritten as 



6iy 



5Vn 



C (l-27^) + 27^, 



where we have defined the ratio 



7^ 



(22) 



(23) 



The splittings linearly increase with C- This linear depen- 
dence is verified using the numerical frequencies computed 
for model Ml and M2 and is shown in Fig. [HI 

We stress that for model M2 the ratio TZ imposed by 
our choice of rotation profile is quite small and the approx- 
imation 

is in agreement with the numerical splitting values. 

We find that the rotational splittings of p-m modes as 
well as g-m modes are dominated by the central layers for 
model M2. For the most g-dominated modes, C ~ 1 and 
the rotational splittings in that case directly give half the 
mean core rotation. 
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Fig. 10. C, = IcoYo/I as a function {v — vq)/ /S.v for alH = 1 
modes for model Ml (top panel) and for model M2 (bottom 
panel), is the frequency of the closest radial mode for 
the 1=1 mode with frequency v. For model Ml, blue dots 
correspond to modes with 9 < v/ /S.v < 15; magenta dots 
to modes with v/Av > 15 and green dots to modes with 
V I Av < 9). For model M2, blue dots correspond to modes 
with 8 < V I A.V < 12; magenta dots to modes with v j Av > 
12 and green dots to modes with v j Av < 8). 



This is not true for less-evolved models such as model 
Ml for which the approximation Eq. ([24| is not sufficient 
(Fig. [S]) . The contribution from the envelope layers is not 
negligible for the rotation profiles we have assumed and the 
surface rotation must be taken into account as in Eq. (|22p . 
This is explained by the fact that the p cavity extends quite 
deep downward inside the model where the rotation already 
shows a sharp gradient toward the center. 

The interest of Eq. (1^^ is that C as well as Si/ and Si/m ax 
can be obtained from observations. iMosser et al.l (I2012b llah 
have shown that it is possible to identify the m = modes 
and their nature (p- or g-dominatcd mixed). We show here 
that it is also possible to attribute them a C value from 
the observations. We use as proxy for the measure of the g 
nature of the mixed mode the distance of the £ — 1 mode 
frequency from that of its companion radial mode: v — vq. 

Figure [TU] shows the ( values ioi £ — 1 modes as a func- 
tion of — ^0 for models Ml and M2 where i^q is the fre- 
quency of the radial mode closest to v. The C values all 
follow the same curve. Due to this monotonic dependence, 
it is possible to estimate the value of C from observations 
{i.e., from the knowledge oi v — i^q). The dispersion of the 
( values about the mean curve in Fig. [TU] introduces uncer- 
tainties on the induced C value. This nevertheless offers a 
mean to derive not only 5fmax ~ hence the mean core rota- 
tion - but also the ratio TZ, hence a measure of the mean 
rotation gradient from observations. 



6. Physical interpretation of the observed splittings 

In this section, we show that rotational splittings can be ex- 
pressed as a function of observables. Indeed, this will further 
permit us to derive physical properties related to rotation 
by using observations only (see Sect. [7]). 

We interpret physically the two components of the ro- 
tational splitting (fi)core and (r2)cnv as the rotation angular 
velocities averaged over the time spent in the g and p cav- 
ities, respectively. We checked numerically with our stellar 
models and numerical frequencies and eigenfunctions that 
(r2)coro and (r2)onv are well approximated by mean rotations 
defined as 

n^-f n{x) dr. (25) 



The time spent by the mode is approximately given by 



fcr da; , 



(26) 



()Unno et al.lll989l and references therein) where a is the 
normalized frequency of the mode. The radial wave number 
of the wave normalized to the stellar radius, fcr, is given, in 
a local asymptotic analysis, by 



1 1 



S f 



(27) 



(Osaki"1975'; lUnno et al.l |1989t) where TV^ and 5"^ are the 
normalized squared Brunt- Vaisala and Lamb frequencies, 
respectively. As will be seen below, the level of approxima- 
tion is sufficient for our purpose. 

In the g cavity, ct^ <^ Sf , therefore 

^ {N' - (28) 
ax 

then the time spent in the resonant g cavity becomes 



(iV2 - ^2)1/2 



dx 



(29) 



Note that this provides the number of nodes in the g cavity 
as Ug that is defined as 2Tmg/a = Tg. 

Similarly, in the acoustic cavity in the envelope where 
CT^ » A^^, Sf, one can write 



Tp = 2 



dx 



(30) 



Note that the time Tp and Tg are defined in units of the 
dynamical time idyn = (G'M/i?^)- 1/2, 

The expressions for the mean core (resp. envelope) ro- 
tation become 



If N 

^core ^ / ^(^X^ dx , 

'^9 ^corc 



ciiv — / fl(x) 



dec 

Cs 



(31) 
(32) 



The theoretical values of fionv and ricoro, Eqs. pip and 
, are computed for both equilibrium models Ml and 
M2. We find 17corc/27r = 0.88 fiRz, Cl^nv/^TT = 0.157 /.tHz 
then 5z/i„ax,th = 0.44 fiRz, TZ = 0.178 for model Ml and 

i^coro/27r — 0.856 /iHz, i7onv/27r = 0.042 /xHz, 5fmax,th = 
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Table 4. Typical values for the quantity 4ao fEa. ([M| ') 
for stars on the RGB and in the clump. Av and i^max are 
given in /iHz and AH in seconds. The values in parenthesis 
corresponds to y„o evaluated at i^max i 3Az/. 
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Fig. 11. Top: C as a function of v j Av. Black open sym- 
bols: numerical results (Eq fTU|) . Red dots: C approximated 
by (Eql2Sl) with x = 2-5 y cos(7r/(Anz^)) and y = v/Av. 
Bottom: Su/Su-a^sj^ as a function of v j IS.v. Black dots: nu- 
merical values (Eg fTTj) . Red dots: approximate expression 
for Jiz/Ji'niax (Eqllll) with 7^ = 0.1785 and C given by its 
numerical values. Blue crosses: bv jbv^^y^ given by Eg 1371 
with 7^ = 0.1785 and x given by Eg [51 

0.428 ^Hz, 7^ = 0.05 for model M2 where the subscript th 
stands for theoretical. 

We then computed the rotational splittings from the 
approximate expression Eq. (j22p using these figures with C, 
computed with the associated eigenfunctions (Eg flOl) . As a 
result, the theoretical dependence given by Eq. (|22l) per- 
fectly matches the numerical curve bv jbv^^-^ and f2corc ^ 
(^^)corG, f^cnv ~ (f^)env This cau bc sccu iu Fig. [11] for 
model Ml where one compares the behavior of the numer- 
ical splittings (5i^/(5i^max with v j Av to the splitting values 
obtained by using Eq. ((22)) . 



One can go to a further level of approximation by de- 
riving an approximate expression for the mode inertia. 
Appendix IA.3I shows that one can derive an approximate 
expression for C, as 



1 



1 + ao 



where 



"0 V 



(33) 



(34) 



For convenience we have defined y — v j Av and a dimen- 
sionless constant as 

ao = Ai/An, (35) 

with An is the period spacing for g-m modes 



An = 



27r^ 



A 



(^2 _ ^2)1/2 



da; 



(36) 



For models Ml and M2, we respectively have ao = 2.2 10 ^ 
and ao = 6.0 lO^*. 

The behavior of the ratio C as a function of v/A.v is 
computed with Eq. ([55)) and Eq. (|M)) and is compared to 
that of the numerical one in Fig. [TT] (top) . The cosine term 
gives rise to the oscillation of C, while the decrease of the 
minimum values of C with frequency is driven by the ratio 
Tp/Tg. ^ is computed using Eq. (|55)) and Eq. ([M)) where 
the factor 2 is replaced by a factor 2.5 which fits better the 
numerical results. This difference is explained by the fact 
that we took the factor / (Appendix IA.3[) equal to unity 
while it is in reality smaller. For model Ml the difference 
between the large separation /S.v and the equivalent quan- 
tity computed over the p resonant cavity is of the order of 
/ = 0.8. This increases the constant in x from 2 to 2.24, 
closer to the the adopted 2.5 value. As mentionned also in 
Appendix lA. 31 the amplitude ratio Eq. (IA.20p is further in- 
creased if one takes into account the effect of the evanescent 
zone between the outer p and inner g resonant cavities of 
modes. As a consequence, this leads to an increase of the x 
term contribution. Actually, some departure from asymp- 
totic is expected for the p part of red giants modes. It is 
then striking that the numerical results agree that well with 
Eq. ([33]). 

An expression for the normalized rotational splitting as 
a function of the observable y = v j Av is then obtained by 
combining Eq. ([22]) and Eq. ([33]) so that 



bv 



5uu 



1 - 2n 

1 + "0 



27^, 



(37) 



with X given by Eq. p4p. 

The rotational splittings computed with Eq. (|M)) and 
Eq. p7p are compared with the numerical values for model 
Ml in Fig. [11] The behavior of both theoretical and numer- 
ical quantities agree quite well. An illustration is provided 
in Se ct. [7]with the red giant star KIC 5356201 (jBeck et all 
'2012) for which such a procedure has been used. 

The ratio Sv/SvYnax as a function of {ly — 1^0) /Sv 
is nearly independent of rotation, in particular when 
7?. <C 1. Maxima of the C oscillation are obtained for 
cos(7r/(ao y)) = 0, then Cmax — 1- Minima are defined for 
cos(7r/(ao2/no)) = ±1, «-e., ao2/no = 1/^*0 (for any integer 
no) then 



1 



"0 Vn 



(38) 



Accordingly, the contrast rj defined as the ratio of mini- 
mum to maximum splittings is obtained from Eq. ([37]) using 
Eq. ([38]) and Eq. ([X23| as 



1 - 27^ 
1 + 4 ao yl 



2TZ. 



(39) 
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Fig. 12. Top; Observed rotational splittings as a func- 
tion of vlb.v for the Kepler star KIC 5356201 (Beck 
et al. 2012). For this star, the mean large separation is 
Aj/ = 15.92/1 Hz while the mean g -mode spaging is 
An = 86.11s and the maximum rotational splittings is 
— 0.405nHz. Bottom: Rotational splittings normal- 
ized to ^fc'max are folded with (y — vq)/ /S.v. 



As 4ao i/ng decreases with increasing integer no, the con- 
trast 5fmin/<5i^max iucrcases with frequency for a given ro- 
tation profile. This can be seen quite clearly for model M2 
in Fig. [5] for instance. 

Typical values of the quantity 4ao 2/n ^"^^ different types 
of observed red giants are listed in Ta b. ^ according to the 
analysis done bv lMosser et al.l (|2012b[ ). For simplicity, yno is 
evaluated at v = i^max- The first two lines are given for red 
giants corresponding well to models Ml and M2 respec- 
tively discussed in the previous sections. More generally, 
this quantity decreases from about 4 at the bottom of the 
RGB down to 0.06 for red giants stars located on the high- 
est part of the RGB for which stochastically oscillations are 
detected. It amounts to roughly 0.4 for the clump stars. 
Note that for the most evolved red giant stars, 4ao is 
quite small. In that case, the contrast rj varies as 



5Vn 



1 - 4 ao yno (1 - 27^) . 



(40) 



and the influence of the rotation on this contrast through- 
out TZ can become negligible. 



7. The illustrative case of the red giant KIC 
5356201 

Differential rotation of the red gi ant KIC 53 56201 observed 
by Kepler has been analyzed bv lBeck et al.l (2012). As the 
inclination of the star has an intermediate value, the three 
components of the dipole mixed-mode multiplets are clearly 
visible, that allows a precise determination of the rotational 
splittings. With z^^ax = 209.7 ± 0.7/iHz, Ai^ = 15.92 ± 
0.02 /iHz and a period spacing AH — 86.2 ± 0.03 s, deter- 
mined according to the method presented in iMosser et al.l 
(|2012b[ ). the star lies at the bottom of the RGB, with a seis- 
mically inferred radius of about 4.5 R(?). Th e minimum rota- 
tional splitting given bv lBeck et al] l|2012[ ) is 0.154 ±0.003 
/iHz; the authors also mention that the average contrast be- 
tween the maximum and minimum v alue of the spli t tings is 
1.7; the maximum value inferred bv iMosser "eTall (l2012bD 
is 0.405±0.010 /tHz. With such global seismic parameters, 
one infers that this star is slightly more evolved than model 
Ml but can be approximately represented by this model. 

The observational rotational splittings are displayed as 
a function of v / (y is the frequency of the centroid mode 
of each identified multiplet) in Fig. [T^] (top). The folded 
splittings are shown in Fig. [T^] (bottom). The global be- 
havior of the curve shows the same characteristics as the 
theoretical equivalent in Fig. [S] for model Ml. Only the 
value of the minimum at {v — vq)/ lS.v is smaller than 0.4 
(Fig. [12]) for the observations compared with 0.59 for model 
Ml (Fig in]). The difference is due to different values of the 
parameters ag (ao = 1.37 10"'^ for the Kepler star). 

Using the observed values for i^max, Ai/ and AH for 
the star KIC 5356201, we find 77 = 0.49 (Eq. 
when evaluated for the mode with frequency v//S.v — 
13.7(~ z^max/Ai^). This is close to the observed contrast 
5v/5vmax = 0.43 evaluated for the same mode. Note that if 
one corrects 77 with the factor / mentioned above, taking a 
typical / = 0.8 as for model Ml, one obtains r] = 0.45 even 
closer to the observed value. 

In Fig. [12] (top), the observed rotational splittings nor- 
malized to their maximum value (Jfc'max are plotted as a 
function of C obtained using Eqs. (j33]) and ([M]) . The linear 
dependence of the observed rotational splittings with C, is 
clear. One single parameter, 27?,, enters the expression for 
5v/5vmax (Eql21]) and can then be estimated to range be- 
tween 0. and 0.1 from Fig. [T3] (top). This indicates that for 
this star the ratio between the average core rotation and 
that of the convective envelope is larger than 20. 

We next compare the theoretical rotational splittings 
computed using Eqs. ([M|) and (j37p with the observed ro- 
tational splittings in Fig. [13] (bottom). The agreement is 
quite acceptable and validates our theoretical derivations 
of Sect, [a 



8. Conclusion 

We have investigated the properties of the rotational split- 
tings of dipole {£ = 1) modes of red giant models. The 
rotational splittings are computed with a first-order pertur- 
bation method because we focused on slowly-rotating red 
giant stars. We considered two models in two different H 
shell-burning evolutionary stages. We find that such modes 
are either g-m modes (inertia and rotational splittings are 
fully dominated by the core properties) or p-m modes (in- 
ertia is almost equally shared between the core and the en- 
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Fig. 13. Top; Observed rotational splittings as a function of 
C computed with Eql33] and EqjMl for the Kepler star KIC 
5358201 (blue dots). The sohd curves are C(l - 27^) + 27^ 
with 7^ = (black) and 27e = 0.1 (red). Bottom: Observed 
rotational splittings (blue filled dots)as a function of v//S.v. 
Theoretical rotational splittings computed according to 
Eg 1571 and EqlM] (black open dots connected with a dashed 
line) with 7^ = 0. 



velope). For g-m modes, the rotational splittings are fully 
dominated by the average core rotation. As the time spent 
by all considered g-m modes in the core is roughly the same, 
the resulting mean core rotation corresponds to the rota- 
tion in the central layers weighted by the time spent by the 
modes in the core. For a radiative core, if the rotation is 
decreasing sharply from the center, the mean core rotation 
is a lower limit of the central rotation. For p-m modes, the 
rotational splittings remain dominated by the core rotation 
but the contribution of the envelope is no longer negligible. 

We find that the rotational splittings normalized to their 
maximum value linearly depend on the core contribution of 
mode inertia. The slope provides the ratio of the average 
envelope rotation to the core one {TV). Mode inertia is uni- 
vocally related to the frequency difference between observed 
£ — 1 modes and the frequency of the closest radial mode. 
As this last quantity is measurable, one can use the simple 
relation between the core inertia and this frequency differ- 
ence to determine the core inertia. This knowledge together 
with a measure of the rotational splitting then provides a 
measure of TZ. 

As a step further in the modeling, we used asymptotic 
properties to show that the core contribution to mode in- 
ertia relative to the total mode inertia depends on the 
frequency normalized to the large separation through a 
Lorentzian. This provides a theoretical support for the use 
of a Lorentzian profile for me asurin g the o bserved split- 
tings for red giant stars (iMosser et al. 2012al ). This also led 
us to find that the behavior of the rotational splittings of 
t = \ modes with the frequency for slowly rotating core 



red giant stars depends on three parameters only. One is 
the large separation Ai/, the second is the g-mode period 
spacing AH, both can be determined observationally and 
can therefore provide again a measure of the third one, TZ. 

The contrast between the minimum {dvrain) and max- 
imum (Jt'max) values of the splittings can be evaluated at 
the frequency of maximum power. It depends only on the 
above 3 parameters. For a core rotating much faster than 
the envelope (negligible TV), we obtain a relation between 
('^t'min), ((^t-max), and AH. This can provide ((J^max) 
(hence the average core rotation) in case when the g-m 
modes have too small amplitudes to be detected and only 
('^I'min) is available. 

Based on the above theoretical developments, we find 
that the Kepler red giant star KIC 53656201 observed by 
Kepler (Beck et al. 2012) rotates with a rotation in the 
convective envelope Vlcz < 0.05 (fi)corc < ^{r = 0) so 
that the core is rotating more than 20 times faster than the 
envelope. 

The present understanding of the properties of rota- 
tional splittings is a contribution to the effort to decipher 
the effect of slow rotation on red giant star frequency spec- 
tra and establish seismic diagnostics on rotation and trans- 
port of angular momentum. For more rapidly rotating red 
giant stars, non-perturbative methods must be used. This 
will be the subject of the third paper of this series. 
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Appendix A: Properties of eigenfunctions and 
mode inertia 

A.l. Displacement eigenfunctions 

The displacement eigenfunctions computed for two modes 
of the selected pattern (see Tabled]) for model Ml are shown 
as a function of the normalized radius r/ R in Fig. lA.ll They 
correspond to the modes with the largest and smallest max- 
imum amplitudes in the pattern. The inner part is domi- 
nated by the horizontal displacement Z2 and oscillates with 
a large number of nodes, typical of a high-order gravity- 
mode. The largest maximum amplitude corresponds to the 
most g-dominated mode whereas the smallest maximum 
amplitudes arise for the p-m modes vi and vq. 

The maximum amplitude of Z2 occurs deep in the g- 
cavity, at the same radius for all modes of the pattern. The 
region of non-negligible amplitude defines the radius of a 
seismic rotating core which is found here independent of 
the mode {r/R ~ 0.02) and far smaller than the upper 
turning radius of the inner gravity resonant cavity {x2 
0.08, Table [3|) 



N 




A. 2. Behavior of /3 and Pcorc with C 

For i — 1 modes of red giants, the term Z2^i in 13 (Eq. (fT3|) ) 
plays almost no role because Z2Z1 ^ z| in the core and 
Z2Z1 <^ zf in the envelope (see Fig. lA.ip . As a result, we 
have 



1 

.1-11 
2 / 



2 I 2 _ 1 
^1+^2 2 ^' 



dx 



dx 



(A.l) 

(A.2) 
(A.3) 



where ^ is defined in EqlTOl The linear dependence of /3 
with ^ is verified in Fig. IA.21 Furthermore, for all modes, 
z| ^ z^, 2ziZ2 in the g- cavity (see Fig lA.ll) then /3coro,ni ~ 



Fig. A.l. Displacement eigenfunctions as a function of the 
normalized radius r/R, in blue the horizontal component 
to the displacement eigenfunction Z2 (Eq. [T]) and in red 
the radial component to the displacement eigenfunction 
zi (Eq. |6j). From top to bottom i — 1 modes 1 and 3 of 
Table, d 



/3coro where for I = 1 modes, we derive 
1 



hence 



/3c 



/ 
1 

21 



7^ + r2 _ 1 :rA — 



dx 



^2 2 2 

1 Ir. 



X 

1 



2 / 



1-C- 



(A.4) 
(A.5) 

(A.6) 
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Fig. A. 2. Top: Pcnv and /3core as a function of i^/Az^ for 
model Ml. Bottom Same as top for the ratio /3corc as a 
function of C (black open dots). The approximation (3 = 
1 — (1/2) ( using the numerical values of C, is represented 
with red crosses. 



Numerical values for model Ml confirm that /3corc in- 
creases linearly with C with a slope 1/2 fFig. IA.2"|) . For g-m 
modes (C ^ 1), /i^corc dominates with a nearly constant value 
of 0.5. P-m modes correspond to the teeth of the saw type 
variation of P^nv and the ratio /3onv//3coro 0.25 (Fig. lA.2[ ). 



A. 3. An approximate expression for 

This section determines an approximate expression of ^ = 
Icore/I as a function of v/ Aiy. The derivation is based on 



resul ts of an asymptotic method developed by iShibahashi 



(Il979i) to which we refer for details (see also lUnno et al 
1198' 



— The envelope ('-^p-propagative cavity) is characterized 
by z| <^ zf. Using Eq. (16.47) from Unno et al., it 
is straightforward to derive the following approximate 
expression 



2 da; / 



1 — sin (TT dr 



(A.7) 



The constant c can be determined by the condition 
1 at the surface. 



with 



2 r 

— / fcr dx' . 



1 



r2 - 5? 



1/2 



(A.8) 



(A.9) 



where we have assumed cr^ >> iV^ in Eql^Tl Note that 
in the process of deriving the amplitude of arising 



in front of the sinusoidal term in Ea lA.7[ one can ne- 
glect Sf in front of in the expression for fcr (i.e. 
fcr ojcs). However this is not the case when fcr is 
in the phase of the sinusoidal term where we keep the 
expression Eq lA.9l 

The inertia in the envelope can then be approximated 
as: 



cnv 
2 



dx 



X 



C ( cos (cTTp) — 1 

— Tp 1 

<J \ (TT, 



P 



ip , 



(AlO) 
(All) 



where we have defined 



I X3 

and the mean large separation is 

10 



Aiy= {2 



(A13) 



The factor / is of order unity and represents the differ- 
ence between the integration from xs and from the cen- 
ter. We take / = 1 unless specified otherwise. The last 
equality in Eq. (|A.lip is obtained assuming aTp ^ 1. 

The core (~g-propagative cavity) is characterized by 
^ zf. Again, the asymptotic results lead to the fol- 
lowing expression 



dx a" / 
X a \ 



1 — sin (TT 1 dr 



(A.14) 



where a is a constant that is determined by the resonant 
frequency condition between the p and g cavities and 



t{xi,x) 



2 



kr dx' . 



and we have used 



ax 



2S1/2 



(A15) 



(A.16) 



Hence recalling that aTg 3> 1 , the inertia in the core can 
be approximated as 



^corc ^ 1 1 

(7 



cos (cTTg) — 1 



(A17) 



where we have defined 



r(xl,x2) = ^ / 



^2)1/2 ^^^^^ 
X 



— The ratio q = Icnv/Icorc is then approximated by 



\aJ Tn 



(A.19) 



— We obtain the ratio c/a (from Unno et al's Eq. (16.49) 
and Eq. (16.50)) as 



c 
a 



C0s(crTg/2) 
COs((TTp/2) 



±2 cosf^ 
V 2 



(A.20) 
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where we have used the fact that (jTp ^ 2rip7r. A ex- 
ponential term is present in Unno et al's expression 
with the argument being an integral over the evanes- 
cent region between the p- and g-cavities. As this re- 
gion is quite narrow in our models for the considered 
modes, the exponential is taken to be 1. Nevertheless 
the width of the evanescent region depends on the con- 
sidered mode and in some cases, for accurate quantita- 
tive results, it might be necessary to include effects of 
the evanescent zone with a finite width. 
The ratio q is then eventually approximated by 



4 cos^ r-^ ^. (A.21) 



For the relative core inertia ^ = Icoic/I, 

(1 + 4 c»^(fi)i)-', (A^22) 

We now use the approximate expressions Eq. (jA.12|) and 
Eq. (jA.18|) in order to derive for the ratio Tp/Tg in terms 
of observable quantities 

^«laoy^ (A.23) 

^9 / 

where for convenience we hase defined y — v/Av and 
ao = Ai^ An, (A.24) 
with the period spacing for g modes 



We also write 



so that we obtain 



= ^ , (A.26) 



1 + ao X // 1 + ao X 



X 



2 4- cos f-r^) . (A.28) 
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